Comment on "Soliton ratchets induced by excitation of internal modes" 
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Very recently Willis et al. [Phys. Rev. E 69, 056612 (2004)] have used a collective variable 
theory to explain the appearance of a nonzero energy current in an ac driven, damped sine-Gordon 
equation. In this comment, we prove rigorously that the time-averaged energy current in an ac 
driven nonlinear Klein-Gordon system is strictly zero. 
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Recently several papers have been published trying 
understand soliton ratchets (see, for example, Refs. 
EJ [J S El and for a recent review chapter 9 in Ref . 
pp. 343-364). This phenomenon is a generalization of 
the ratchet effect p| to spatially extended systems, and 
manifests as a unidirectional motion of a soliton induced 
by zero-average forces. A paradigmatic example is the 
driven, damped nonlinear Klein-Gordon equation: 

4>,u(x, t) - 4>, xx (x, t) = -U' [<(>(x, t)} + f(t) - 0<f>, t (x, t), 

(1) 

where g tZ — dg/dz, f(t) is a periodic field with period 

T and zero time-average [i.e., 1/T J Q T dtf(t) = 0], /? > 
is the dissipation parameter determining the inverse 
relaxation time in the system, and U'(z) is the derivative 
with respect to z of the potential U(z). In this comment, 
we will assume that the potential U(z) is periodic with 
period A, and presents minima at Zj = zq + jX, with 
j G Z. The ac driven, damped sine-Gordon equation 
considered in Ref. [6j is a particular case of this more 
general problem, with U (z) = 1 — cos(z) and 



f(t) = - [ei cos(wt) + e 2 cos(2ujt + 0)} 



(2) 



To fully specify the mathematical problem, the partial 
differential equation must be amended by both ini- 
tial conditions for 4>{x, 0) and <f>,t(x, 0), and boundary 
conditions for lim x ^± 00 4>{x, t). Several boundary condi- 
tions can be imposed to have a well-posed boundary value 
problem. For instance, in the absence of the periodic field 
f(t), it is possible to choose the fixed boundary condi- 
tions: lim x _ >+00 4>(x, t) = zi and lim^-oo <f>(x, t) = z m . 
In the presence of f(t), the fixed boundary conditions 
become incompatible with Eq. QJ, and they are usually 



replaced by the aperiodic boundary conditions: 

lim (f)(x,t) = lim <f){x,t) + XQ, (3) 

x — >+oo x — > — 00 



lim cj) x (x,t) 

X — ' + 00 



lim (f> tX (x,t), 



(4) 



where Q G Z is the so-called topological charge. The 
discrete version of these periodic boundary conditions are 
also the most used in the numerical solution of Eq. (JTJ 
(see, for example, Refs. and [j|). 

It can be derived from the continuity equation that 
the energy current density generated by the field <j)(x, t) 
in the absence of damping and external forcing is given by 
j(x, t) = —4> t t(x,t)4> tX (x, t) and, consequently, the energy 
current reads 



J(t) = 



dx 4>,x{x,t)(j> : t{x,t). 



(5) 



The time-averaged energy current (J) is defined as the 
limit 



(J) 



lim 



1 



dt J(t). 



(6) 



In Ref. 0, it has been proved by symmetry considera- 
tions that a necessary condition for the appearance of 
a non-vanishing time-averaged energy current is that ei- 
ther the potential presents broken spatial symmetry, or 
the field f(t) violates the symmetry property 



/(*+!) =-/<"■ 



(7) 
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or both simultaneously. Following this idea, a collec- 
tive variable approach has been developed in Ref. @ 
for the ac driven, damped sine-Gordon equation with a 
field of the form © which leads to a non- vanishing time- 
averaged energy current. The purpose of this comment is 
to prove that the time-averaged energy current, (J), of a 



2 



driven, damped nonlinear Klein-Gordon equation of the 
form Q is necessarily zero and, consequently, the above 
mentioned results in Ref. Q must be erroneous. 

To prove that (J) = 0, firstly we will obtain an ordi- 
nary differential equation for the energy current J it). In 
order to do that, we differentiate with respect to time 
Eq. O, resulting 



/+oo 
-oo 

(8) 

By making use of Eq. 0) in the second term on the right- 
hand side of the above expression, it is straightforward 
to write it in the form 



j(t) 



-1 lim [<t). t (x,t)} 2 - lim [<j>. t (x,t)} 2 

Z I a;— >-foc - 



x — X. 



lim [(f> x (x,t)] - lim [<f>. x (x,t)] 2 

Z I x— »+oc x— oc 



+ lim U[4>(x,t)]- lim U[(j>(x,t)] 



x— *-\-oo 



-0J(t) 



x^ — oo 



lim < 

x— »+oc 



(x, t) — lim (^(a:, £) 
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(9) 



Differentiating Eq. with respect to t, it is easy to see 
that the first term between brace brackets on the right- 
hand side of Eq. © is equal to zero. From Eq. (0J, it 
follows that the second term between brace brackets on 
the right-hand side of Eq. © is also equal to zero. The 
two terms of Eq. © containing U [<f)(x, t)] also cancel 
each other due to the boundary condition © and the 
periodicity of U(z). Thus, from Eq. |JS} we finally obtain 



j(t) = -0J(t) - XQf(t). 



(10) 



It is important to emphasize that Eq. I(10|) is a direct con- 
sequence of Eq. ||TJ and the boundary conditions |J3J and 
PJl. Therefore, it is an exact result valid for any periodic 
potential U (z) of the type described in the paragraph be- 
low Eq. QJ, and any external field f(t). Equation ltTU|) 
appears in Ref. |(J as an approximate result obtained af- 
ter neglecting the dressing due to phonons. 
The general solution of Eq. (fT~0f> is 

J(t) = J{0)e~ pt - XQ I alt' e^'"*') f(t'), (11) 



and making use of the definition of the time-averaged 
energy current in Eq. ©, it results 



w = 



dt f(t) 



\q r 



fir J 



dt f{t)e-^ T -^ 



(12) 



The first limit in the above expression is obviously equal 
to zero. The second one is also equal to zero as the 
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FIG. 1: Evolution of the energy current J(t) corresponding to 
the sine-Gordon equation with the biharmonic field in Eq. @ 
and the initial conditions corresponding to a static kink-likc 
solution of the unperturbed sine-Gordon equation [i.e., f(t) = 
/3 = 01. The parameter values are the same as in Fig. 3 of 
Ref. 0: /3 = 0.15, ei = 0.16, e 2 = £i/a/2, lu = 0.25, and 
9 — 1.61 — 7r. The points represents the results obtained by 
solving numerically the sine-Gordon equation and the solid 
line the analytical expression in Eq. 1)14^ . The agreement is 
obviously excellent since Eq. (1141 is exact. The time-averaged 
energy current (J) is zero. 



external field is periodic with zero time-average. The 
integral appearing in Eq. 1)12(1 can be bounded using the 
fact that 



dt /(t)e-« r - f) 



< 



< 



dt \f{t)\e-^ T - t) 



T [ 



,-0r\ 



(13) 



where f m is an upper bound of \f(t) \ and, thus, the third 
limit in Eq. I|12|l is also equal to zero. We conclude that 
(J) = for any periodic potential U(z) of the type de- 
scribed in the paragraph below Eq. Q), and any bounded, 
zero time-averaged periodic field f(t). 

To corroborate this result, we have solved numerically 
the sine-Gordon equation for the case of the biharmonic 
field in Eq. We have considered the initial conditions 
for <p(x, 0) and (p,t{x, 0) corresponding to a static kink-like 
solution of the unperturbed sine-Gordon equation [i.e., 
f(t) =0 = Q], so that, X = 2tt, Q = 1, and J(0) = 0. 
For this particular choice, the energy current in Eq. I|ll|) 
reads 

27rei 



J(t) 



+ LO Z 

2ire 2 



-0e 



-I3t 



- 0cos(ujt) + ujs 

0t {0 cos6> + 2w sinf 
+ 2u sin(2wi - 



i(wt)] 



2 + Auj 2 1 

+0cos(2ujt + 9) +2wsin(2u>i + 0)]. (14) 

Figure n shows the perfect agreement between the above 
expression and the numerical simulation of the sine- 
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Gordon equation. This agreement is not surprising at 
all since Eq. (|14f> is exact. The time-averaged energy 
current (J) is zero. 

It is important to emphasize that the result in this 
comment is not applicable when the external field not 
only depends on t but also on x. In that case Eq. 
cannot be obtained and, in principle, it is possible to 
observe a non- vanishing time-averaged energy current. 
A field of this kind has been considered in Ref. p) , where 
f(x,t) = E(t) + £(x,t), with E(i) being an ac field with 
zero mean and i) a Gaussian white noise. 
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